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hancement. The FCNC decays B° S (B°) -> fi+fj,' [1( 
occur in the SM only through higher order diagrams. 
The SM expectations for these branching fractions are 
B(B° S -> fi + n~) = (3.42 ± 0.54) x 10" 9 and B(B° -> 
/x + M~) = (1.00 ± 0.14) x 10" 10 0, which arc one order 
of magnitude smaller than current experimental sensitiv- 
ity. Previous bounds, based on 1.3 fb _1 and 364 pb _1 
arc B{B° S -> ^ + pr) < 1.2 x 10~ 7 and B(B° -> < 
5.1 x 10~ 8 at the 95% C.L., respectively di]. 

Enhancements to B®(B°) — > n + n~ occur in many new- 
physics models. In supersymmetry (SUSY) models, con- 
tributions from diagrams including supcrsymmctric par- 
ticles can increase B(B°(B°) — > ,u + /i~) by several orders 
of magnitude at large tan/3, the ratio of vacuum expec- 
tation values of the Higgs doublets [B|. In the minimal 
supcrsymmetric standard model (MSSM), the enhance- 
ment is proportional to tan 6 /3. Global analyses includ- 
ing all existing experimental constraints suggest that the 
large tan/3 region is of interest In contrast, SUSY 

R-parity violating models [|| and non-minimal flavor vi- 
olating models [9(] can both enhance B® — ► /J + /i~ and 
B° — > ^t + /i _ separately even at low tan/3. In the absence 
of an observation, limits on B(B° — ► ^t + /i~) are comple- 
mentary to those provided by other experimental mea- 
surements, and together would significantly constrain the 
allowed supersymmetric parameter space. For example, 
if the lightest neutralino in SUSY models is a cold dark 
matter (CDM) particle, B(B® — > fi + fi~) and constraints 
on the amount of CDM in the universe from cosmic mi- 
crowave anisotropy measurements can be exploited in 
this way @, 0, @l- Then, for instance, in minimal super- 
gravity (mSUGRA) models limits on B{B° S — > n + will 
correspond to bounds on supcrpartner particle masses 
that are beyond the sensitivity of the corresponding di- 
rect searches for those particles in colliding-beam exper- 
iments [H . In general, the search for these rare decays is 
central to exploring a large class of new-physics models. 

This measurement uses 2fb _1 of integrated luminosity 
collected by the upgraded Collider Detector at Fermilab 
(CDF II) and supersedes our previous measurement using 
364 pb -1 0]. The sensitivity of the analysis is improved 
significantly by increasing the integrated luminosity of 
the event sample, using an enhanced muon selection, em- 
ploying a neural network (NN) discriminant to separate 
signal from background, and performing the search in a 
two dimensional grid in dimuon mass and NN space. 

A detailed description of the CDF II detector can be 
found in Ref . [l(| • Charged particle tracking is provided 
by a silicon microstrip detector surrounded by an open- 
cell wire drift chamber immersed in a 1.4 T solenoidal 
magnetic field. This system provides precise vertex de- 
termination and momentum measurements for charged 
particles in a pseudorapidity range \r/\ < 1.0, where 
j] = — ln(tan |) and 9 is the angle of the track measured 
with respect to the proton beam direction. Additionally, 
the drift chamber measures the ionization per unit path 
length, dE/dx, for particle identification. Surrounding 



the tracking detectors are electromagnetic and hadronic 
sampling calorimeters arranged in a projective geometry. 
Drift chambers referred to as CMU and CMX are located 
behind the calorimeters to detect muons within \rj\ < 0.6 
and 0.6 < \r/\ < 1.0, respectively. 

Events are recorded for subsequent analysis if they 
have either of two topologies that satisfy the require- 
ments of the online trigger system: CMU-CMU, which 
includes events where both muon candidates are trig- 
gered using the central muon detectors, and CMU-CMX, 
which includes events where one of the muons is triggered 
in the central muon detector and one in the higher pseu- 
dorapidity region. The details of the trigger system and 
selection requirements can be found in Refs. [S, Ell- 
in the offline analysis, the trigger selection is refined by 
applying a series of "baseline" requirements. We select 
two oppositely charged muon candidates within a dimuon 
invariant mass window of 4.669 < m w < 5.969 GeV/c 2 
around the B® and B° masses. Backgrounds from 
hadrons misidentified as muons are suppressed by select- 
ing muon candidates using a likelihood function. This 
function tests the consistency of electromagnetic and 
hadronic energy with that expected for a minimum ion- 
izing particle and the differences between extrapolated 
track trajectories and muon system hits (l2l |. In addi- 
tion, backgrounds from kaons that penetrate through the 
calorimeter to the muon system or decay in flight outside 
the drift chamber are further suppressed by a loose se- 
lection based on dE/dx fl3j . The inputs to the muon 
likelihood and the dE/dx performance are calibrated us- 
ing samples of J/ij) — > D° — > K~tt + and A — > pir~ 
decays. To reduce combinatorial backgrounds the muon 
candidates are required to have transverse momentum 
relative to the beam direction pr > 2.0(2.2) GcV/c for 
CMU(CMX), and \p^\ > 4 GcV/c, where p^ is the 
transverse component of the sum of the muon momen- 
tum vectors. The remaining pairs of muon tracks are fit 
under the constraint that they come from the same three- 
dimensional (3D) space point. To achieve further sepa- 
ration of signal from background, we employ additional 
discriminating variables. As in our previous analysis [J 
these include the measured proper decay time, A; the 
proper decay time divided by the estimated uncertainty, 
\/cr\; the 3D opening angle between vectors p^ and 
the displacement vector between the primary vertex and 
the dimuon vertex, A0; and the 5-candidate track iso- 
lation, I Q. We require that \/a x > 2, A 8 < 0.7 rad, 
and / > 0.50. There are 30 666 dimuon candidates that 
fulfill the above trigger and baseline reconstruction re- 
quirements. The baseline selection reduces combinato- 
rial backgrounds by a factor of 300 while keeping ap- 
proximately 50% of the signal events that are within the 
acceptance (geometric and kinematic requirements) of 
the trigger. Relative to the previous analysis, the ap- 
plication of the muon likelihood and dE/dx selection is 
approximately 97% efficient for the signal and reduces 
combinatorial backgrounds by 35% and constrains back- 
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FIG. 1: Distributions of vn for simulated B® — * fi + fi signal 
and observed sideband events. 



same efficiencies and optimizations are used. 

For measuring efficiencies, estimating backgrounds, 
and optimizing the analysis, samples of B®(B°) — > 
B + — > J/ipK + , and B — > h + h~ are generated with the 
PYTHIA simulation program [l6[ and a CDF ff detec- 
tor simulation. The B-hadron px spectrum and the I 
distribution of the i?-hadrons are weighted to match dis- 
tributions measured in samples of B + — ► J/rpK + and 
B° s -> J/V></> events [T3, E3- 

We use a relative normalization to determine the B® — > 
branching fraction: 



J_ fu 

e N fs 



B(B + ), (1) 



where N s is the number of B® — ► candidate events. 

We observe N + = 11 387 ± 164 B+ -> J/V>A'+ candi- 
dates. This number has been corrected for background 
using sideband subtraction; the contribution of £?+ — > 
J/%pTT + events is negligible. We use 2?(£? + ) = B(B + — > 
J/V>A+ -> n + ii~K+) = (5.94±0.21) x 10" 5 and the 
ratio of B-hadron production fractions f u /fs = 3.86 ± 
0.59 [3]. The parameter a s (a+) is the acceptance of the 
trigger and e s (e+) is the efficiency of the reconstruction 
requirements for the signal (normalization) mode. The 
reconstruction efficiency includes trigger, track, muon, 
and baseline selection efficiencies. The NN efficiency, 
eat, only applies to the signal mode. The expression for 
B(B° — > /i + fj,~) is derived by replacing with i?° and 
the fragmentation ratio with f u / fd = 1- 

The ratios of acceptances a+/a s are 0.297 ±0.020 and 
0.191 ± 0.013 for the CMU-CMU and CMU-CMX trig- 
ger channels, respectively. These ratios are measured 
using simulated events, and the uncertainties include 
contributions from systematic variations of the model- 
ing of the £?-hadron px distributions, the longitudinal 
beam profile, and from the statistics of the simulated 
event samples. The ratio of reconstruction efficiencies 
is e+/e s = 0.89 ± 0.04. Muon reconstruction efficiencies 
are estimated as a function of muon px using observed 
event samples of inclusive J/ip — > // + n~ decays. System- 
atic uncertainties in the efficiency ratio, e+/e s , largely 
cancel with the exception of the kaon efficiency from 
the B + decay. The uncertainty is dominated by kine- 
gions 4.669 < m w < 5.169 GcV/c 2 and 5.469 < m MAt < mati c differences between inclusive J/ip -> n+pT and 



grounds from two-body hadronic B decays to a level 
that has negligible impact on the analysis. A sample of 
B + — > J/ipK + events is collected to serve as a normal- 
ization mode using the same baseline requirements, but 
including a requirement of pr > 1 GeV/c for the kaon 
candidate and constructing the B + — > J/tpK + — > 
vertex using only the muon candidate tracks. 

For the final event selection we use the following dis- 
criminating variables: m aa , A, \/<j\, AG, /, |, and 
the pr of the lower momentum muon candidate. To en- 
hance signal and background separation we construct a 
NN discriminant, v?j, based on all the discriminating 
variables except m w , which is used to define signal and 
sideband background regions. The NN is trained using 
background events sampled from the sideband regions 
and signal events generated with a simulation described 
below. The vn distributions of B® signal and sideband 
background events are shown in Fig. [1] Only part of the 
total number of background and simulated signal events 
are used in order to have unbiased samples to test the 
background discrimination and signal efficiency. 

For the final selection, we define search regions of 
5.310 < m w < 5.430 GeV/c 2 for the B a s and 5.219 < 
< 5.339 GeV/c 2 for the B° around the mass val- 
ues m B o = 5.370 GeV/c 2 and m B o = 5.279 GeV/c 2 [HI, 
respectively. These regions correspond to approximately 
±2.5 times cr m , the estimated two-track invariant mass 
resolution, where a m ps 24 McV/c 2 . The sideband re- 



5.969 GeV/c 2 are used to estimate the combinatorial 
backgrounds in the signal regions. Backgrounds from the 
two-body hadronic B® and B° decays, B — ► h + h~ , where 
h are ^ or K ± , which peak in the B° s and B Q invari- 
ant mass signal region and do not occur in the sidebands, 
are estimated separately. The content of signal regions 
were not unveiled until all selection criteria were final- 
ized. The final selection criteria were determined from 
an a priori optimization, which maximizes sensitivity of 
the expected limit. The kinematics of — > and 
B° — > decays are similar enough that the efficien- 

cies are the same within statistical uncertainties and the 



Bg(B°) — > decays. The efficiency, cn, is estimated 

from the simulation. We assign a relative systematic un- 
certainty on cat of 6% based on comparisons of NN per- 
formance in simulated and observed B + — ► J/ipK + event 
samples and the statistical uncertainty on studies of the 
B° pt and I distributions from observed B® — > J/if-'cj) 
event samples. The NN provides approximately 25% bet- 
ter background rejection for the same signal efficiency 
compared to the previous analysis 0]. 

The expected background is obtained by summing con- 
tributions from the combinatorial continuum and from 
B — > h + h~ decays. The contribution from other heavy- 
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flavor decays is negligible. We estimate the combinatorial 
background by linearly extrapolating from the sideband 
region to the signal region. The B — * h + hr contribu- 
tions are about a factor of ten smaller than the combi- 
natorial background and are estimated using efficiencies 
taken from the simulation, probabilities of misidentify- 
ing hadrons as muons measured in a D a — > irK data 
sample, and normalizations derived from branching frac- 



tions from Refs. [131. 1 1 51] - The two-body invariant mass 
distribution of the simulated B — ► h + h~ candidates is 
calculated from the momentum of the hadrons assum- 
ing the muon mass hypothesis. The background esti- 
mates are cross-checked using three independent control 
samples: /J^fJ^ events, events with A < 0, and 

a misidentified muon-enhanced /i + \i~ sample in which 
we require one muon candidate to fail the muon quality 
requirements. We compare the predicted and observed 
number of events in these samples for a wide range of un 
requirements and observe no significant discrepancies. 

Using an a priori optimization procedure, we found 
that subdividing the signal region into several and 
mass bins to exploit the shape of the mass distribution 
and the higher signal to background ratios for the higher 
vjq values improves the sensitivity by 15% relative to us- 
ing a single bin. The signal region is divided into five 
equal mass bins of 24 MeV/c 2 and three bins delin- 
eated at 0.8, 0.95, 0.995 and 1.0. The backgrounds, effi- 
ciencies, and limits are calculated in each bin separately. 
Summing over the mass bins in each slice of , the cor- 
responding eats are estimated to be 12%, 23%, and 44% 
and the expected SM yields of B® — > events are 

0.08±0.03, 0.15±0.05, and 0.30±0.10, respectively The 
expected yield of B° — > fJ, + fi~ events is ten times smaller. 
Using these optimized selection criteria, we compute an 
expected limit of B(B° S -> < 4.9 x 10~ 8 at 95% 

C.L. The expected limit is a factor of five better than the 
expected limit of the previous analysis Q ■ The limits are 
estimated from Eq. (JTJ) using the confidence level (CLs) 
method of Ref. [15[ to extract the 95% C.L. upper bound 
on N s ; the limit incorporates Gaussian uncertainties on 
the signal acceptance and efficiency as well as the back- 
ground estimates. The number of observed events is com- 
pared with the total predicted background in Table Q] for 
each bin of mass and vn. The uncertainty on the back- 
ground estimate is dominated by the statistical uncer- 
tainty of the sideband sample. The /U + /i~ invariant mass 
distributions for the three different vjy ranges are shown 
in Fig. (5) The observed event rates are consistent with 
SM background expectations. We extract 95% (90%) 
C.L. limits of B{B° S -> < 5.8 x 10~ 8 (4.7 x 10~ 8 ) 

and B(B a -> < 1-8 x lO" 8 (1.5 x 10" 8 ). 

In mSUGRA, branching ratios as low as B(B® — > 
= 5.8 x 10 -8 occur for common gaugino mass 
parameter, m 1 / 2 , below 380 GeV/c 2 at tan/3 = 50 in a 
CDM-allowed co-annihilation region [l7| • In this scenario 
we exclude gluino masses below 925 GcV/c 2 . 

This Letter reports a search for the rare FCNC decays 
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FIG. 2: The /x + fi invariant mass distribution for events sat- 
isfying all selection criteria for the final three ranges of un- 





mass bins 






v N I II 


III 


IV 


V 


A Exp 10.3±0.5 10.1±0.4 9.9±0.4 


9.7±0.4 


9.5±0.4 


Obs 11 9 


10 


9 


5 


B Exp 3.7±0.3 3.7±0.3 


3.6±0.3 


3.5±0.3 


3.5±0.3 


Obs 4 3 


6 


6 


2 


C Exp 0.7±0.1 0.7±0.1 


0.7±0.1 


0.7±0.1 


0.7±0.1 


Obs 1 


1 





1 


B° 








A Exp 11.0±0.6 10.8±0.5 10.7±0.5 10.5±0.5 10.3±0.5 


Obs 15 13 


9 


14 


9 


B Exp 4.0±0.3 3.9±0.3 


3.9±0.3 


3.8±0.3 


3.7±0.3 


Obs 1 1 


5 


2 


4 


C Exp 0.8±0.1 0.8±0.1 


0.8±0.1 


0.8±0.1 


0.8±0.1 


Obs 2 3 


1 









TABLE I: The total number of expected(Exp) and ob- 
served(Obs) background events for the (upper) and B° 
(lower) signal windows. v N bins: A (0.80-0.95), B (0.95- 
0.995), C (0.995-1.0) and five equal sized mass bins, (I-V), 
as described in the text. 



B® — > /-t + yU _ and B° — > /i + /-t _ with 2 fb~ integrated 
luminosity collected in pp collisions at y/s = 1.96 TeV 
using the CDF II detector and employing improved anal- 
ysis techniques. We observe no evidence for new physics 
and set limits that are the most stringent to date, im- 
proving the previous results 0, 0] by a factor of two or 
more. These limits place further constraints on new- 
physics models @, 0, 0, 1, 0, 

and complement direct 
searches for new physics. We expect the analysis sensi- 
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tivity to continue to improve as we include larger data 
sets. 
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